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rj ■ ABSTRACT 

. \ We use a large sample of upper limits and accurate estimates of supermassive black 

■ holes (SMBHs) masses coupled with libraries of host galaxy velocity dispersions, rota- 

tional velocities and photometric parameters extracted from Sloan Digital Sky Survey 
Q , i-band images to establish correlations between the SMBH and host galaxy param- 

( ■ eters. We test whether the mass of the black hole, M,, is fundamentally driven by 

^ \ either local or global galaxy properties. We explore correlations between M, and stel- 

es . lar velocity dispersion cTe, «-band bulge luminosity Li,buigo, bulge mass Mbuigc, bulge 

Sersic index n, bulge mean effective surface brightness (/io, bulge), i-band luminosity 
\ of the galaxy -/ji,gai, galaxy stellar mass M^^gai, maximum circular velocity T4, galaxy 

^ . dynamical and effective masses Mdyn.gai and Mo^gai. We verify the tightness of the 

' M, — CTe relation and find that correlations with other galaxy parameters do not 

VjQ , yield tighter trends. We do not find differences in the M, — CTc relation of barred and 

CSJ ■ unbarred galaxies. The M, — CTc relation of pseudo-bulges is also coarser and has a 

' different slope than that involving classical bulges. The M,— Mbuigo is not as tight as 

, the M, — CTc relation, despite the bulge mass proving to be a better proxy of M, than 

' bulge luminosity, and despite adding the bulge effective radius as an additional fitting 

T— I \ parameter. Contrary to various published reports, we find a rather poor correlation 

T-H . between M, and n (or (/ic, bulge)) suggesting that M, is not related to the bulge light 

^ ' concentration. The correlations between Af, and galaxy luminosity or mass are not a 

1-H , marked improvement over the M, — CTc relation. These scaling relations depend sensi- 

^ ■ tively on the host galaxy morphology: early-type galaxies follow a tighter relation than 

5-H ' late-type galaxies. If Vc is a proxy for the dark matter halo mass, the large scatter of 

the M, ~ Vc relation then suggests that M, is more coupled to the baryonic rather 
than the dark matter. We have tested the need for a third parameter in the M, scaling 
relations, through various linear correlations with bulge and galaxy parameters, only 
to confirm that the fundamental plane of the SMBH is mainly driven by CTe with a 
small tilt due to the effective radius. We provide a compendium of galaxy structural 
properties for most of the SMBH hosts known to date. 

Key words: black holes physics — galaxies: fundamental parameters — galaxies: 
photometry — galaxies: kinematics and dynamics — galaxies: statistics 



1 INTRODUCTION 



The mass M, of supermassive black holes (SMBHs) is 
closely tied to the properties of the spheroidal co mponent 
of ga la xies, such as the bulge l u minos ity, Lbuige JPressler 
1 19891: iKormendv fc Richstonel 1 19951 : iMarconi fc Hunt 
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I2OO3I : iGrahamI [2OO7I : iGiiltekinet all |2009|. hereafter lG09h 
the ste llar velocity d is persion , a jPerrarese fc Merritj 
2OOOI: iGebhardt et all I2OO0I : pKem ainc et al. 200j; 
Ferrarese fc Ford l2005l : |G09|). the mass of the bulee 



i Magorrian et al.l 



Haring fc Rb3 |2004| '1. the central 



light concentration (iGraham et al.l 200 ll ) , the Sersic index 
(Graham fc Driver 2007), the virial mass of the galaxy 
(jFerrarese et al.i i2006i ). the gravitational binding energy 
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l|Aller fc Richstond 1200 if), the kinetic ener gy of random 
motions of the bulge ( Feoli fc Mancinill2009l ). and the stel- 
lar light and mass deficit associated to the cor e ellipticals 
l|Lauer et al.l l2007l : iKormendv fc Benderl |2009|). Most o f 
these relati ons are inter-compared in Novak et alj (|2006| ) 
and in [Gogl . Given the M. — cr relation a nd the correlation 
bet ween a and the circ ular velocity, \4, iFerraresd (|2002l ') 
and lPizzella et al.l l|2005! ) suggested a link between M. and 
Vc (or equ ivalently, with t h e mas s of th e da r k mat ter halo). 
However, ICourteau et all (|2007l ') and [h3 (|2007l 'l pointed 
out that the Vc—u relation actually depends on galaxy 
morphology (or equivalently, on its light concentration) 
thus precluding a simple c onnection between M. and Vc- 
iKormendv fc Bender] l|201ll ) studied a sample of bulgeless 
galajcies and concluded that there is almost no correlation 
between the SMBH and dark matter halo, unless the galaxy 
also contains a bulge. 

Several authors have noted that the residuals of the 
M, — a and M. — L buige relations correlate with the galaxy 
effectiv e radius fe.g.. lMarconi &: Ifuntll200"3 'l. iHopkins et al.l 
l|2007al lbl) suggested the possibility of a linear combina- 
tion between different galaxy properties to reduce the scat- 
ter of the M, scaling laws, heralding the idea of a fun- 
damental plane of SMBHs (BHFP). Many SMBH scaling 
relations could thus be seen as projections of the BHFP 
lAller fc Richstonell2007l : iBarwav fc Kembhavil20"07l ) . A cor- 
relation of M, with more than one galaxy parameter would 
suggest a SMBH growth sensitive to the overall structure of 
the host galaxy. 

The local characterisation and cosmic evolution of 
the M, scaling relations have already been examined 



througii tJieoreti cal models tor tiie coevolution or galajc- 
ies and SMBHs (iGranato et al.ll2004l: IVittorini et al.ll2005l : 



I Hopkins et al. ll2006l : lMonaco et al.ll2007i r These studies have 
revealed that the observed scaling relations could be re- 
produced in models of SMBH growth wi th strong feedback 
from the active galactic nucleus (AGN . ISilk k, Reed 1 199^: 
Cox et aLllioO^ : [Robertson et all l2006al lbl ; iDi Matteo et al.l 
20051 ). In particular, these models predict the existence of 



the BHFP (|Hopkins et all [2007a b. 2009). However, while 
the observed relations can be reproduced by the models, 
thes e still depend on the adopted slope, zero point, and scat- 
ter l|Somervilldl2009l ') which remain ill-constrained. 

In this work, we make use of a large sample of galaxies 
with available M, estimates to improve our understanding 
of the known scaling laws over a wide range of M,, mor- 
phological type and nuclear activity, as well as to test for 
possible correlations of M, with different combinations of 
spheroid and galaxy parameters. 

This paper is organised as follows. The sample of SMBH 
hosts is described in ij2l Their photometric, kinematic, and 
dynamical properties are presented in 33] while the correla- 
tions between M, and the bulge and galaxy properties are 
shown in 34] We discuss our results and conclude in 35] 



2 SAMPLE DESCRIPTION 

The M, values were retrieved from t wo different samples 



the comp ilatio n of M, upper limits bv lBeifiori et al. 
hereafter IB09| ) and the co mpila tion of secure M, by ; 

The M. estimates of IB09] were obtained from Hubble 



il2oog . 
Gog . 



Space Telescope {HST) archival spectra. That sample in- 
cludes nuclear spectra for 105 nearby galaxies with D < 100 
Mpc obtained with the Space Telescope Imaging Spectro- 
graph (STIS) equipped with the G750M grating. The spec- 
tra cover the region of the Ha line and [N ll] AA6548, 6583 
and [S ii] AA6716, 6731 doublets. The nebular-line widths 
were modelled in terms of gas motion in a thin disc of un- 
known orientation but known spatial extent following the 
method of Sarzi et al. ( 2002. ) . ,B09, adopted two different in- 
clinations for the gaseous disc with a nearly face-on disc 
(i — 33°) hosting a larger M, and a nearly edge-on disc 
(81°) harbouring a smaller M,. The two inclinations cor- 
respond to the 68% upper and lower confiden ce lim its for 
randomly oriented discs. We augmented the IBOGI sample 
with the M. upper limits of NGC 2892 and NGC 5921. The 
STIS/G750M spectra for these galaxies were retrieved from 
the HST archive and we have calculated their M, upper 
limits from the nebular line widths following the prescrip- 
tion of IB09] . We include in Sample A the set of 105 g alaxies 
from IBOGI minus 18 galaxies in common with IG09] (15 of 
them with a secure M, and 3 with a M, upper limit derived 
from the dynamical modelling of resolved kinematics). We 
are left with 87 M. upper limits from lBOgI based on nebular- 
line widths. We also include two newly determined M, upper 
limits and the five upper limits deri ved fr om the dynamical 
modelling of resolved kinematics by IGOGI . The resulting 94 
galaxies constitute our Sample A. 

[Gogj collected M, data published up to November 2008 
based on the resolved kinematics of ionised gas, stars, and 
water maser for total sample of 49 galaxies with a secure 
Al, estimate and 18 galaxies with a Al, upper limit. Our 
Sample B is limited to the 49 definite values of Al,. The 
five upper limits derived from the dynamical modelling of 
resolved kinematics already belong to Sample A. The re- 
maining 13 upper limits by G09 arc taken from the nebular 
line width measurements by ISarzi et al.1 (|2002l ). They were 
also measured by IB09] and are therefore already included in 
Sample A. 

All the data relative to Samples A and B ar e list ed in 
Table [T] and [2] respectively. The upper limits by |B09| were 
rescaled assuming Ho = 70 km s~^ Mpc"\ = 0.3, and 
Qa = 0.7 to conform with IGOgI . 

Figure [T] shows the upper limits and the accurate de- 
ter minat ions of M , for the 18 galaxies in common between 
the iBOgI and |G09| samples. The M, estimates by lB09l and 
iGogI are consistent within 1 cr of each other; no systematic 
offset is detected. Thus, nebular line width measurements by 
[b09 (included in Sample A) trace well the nuclear gravita- 
tional potential dominated by the central SMBH, allowing 
a reliable estimate of M, (as those in Sample B). There- 
fore, we shall use the M, upper limits from Sample A to 
study the correlations of M, versus various galaxy parame- 
ters. We adopt for our tests th e cas e of i = 33° which max- 
imises the upper limit on M.. IB09I already compared their 
M, — fje relation w ith those of iFerrarese fc Ford] ()2005] ) and 
iLauer et al.l ()2007] ) to show that their upper limits on M, 
(included in Sample A) are a valuable proxy for the more 
secure determinations of M. (comprised in Sample B). A 
Kolmogorov-Smirnov test on the distributions of M, and cr 
indicates that Samples A and B could be drawn from the 
same parent distribution to better than the 80% confidence 
level. For this reason, we merged Samples A and B into a 
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3.1 Galaxy Photometric Parameters 

We retrieved o and z-band images f rom tlie seventli data 



release (DR7, lAbazaiian et al 
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Figure 1 . Comparison between tlie M, upper limits by IB09I and 
accurate M, me asurements (symbols) and upper limits (leftward 
arrows) bv lG09l and based on the resolved kinematics of gas (filled 
circles), stars (open circles), and water mascrs (open sq uare) . The 
upper and lower edges of the dotted lines correspond to iB09l 's M, 
values estimated assuming an inclination of i = 33° and 81° for 
the unresolved Keplerian disc, respectively. 



joint A+B data set. The combination of Samples A and B 
yields a total sample of 143 M, determinations. The union 
of those two samples proves most valuable for the proper 
statistical assessment of scaling relations with M,. 



3 GALAXY PROPERTIES 

The aim of this study is to determine the strength of the 
correlations, if any, between M, and the properties of their 
host galaxies. For the latter, we have compiled as large a col- 
lection as possible of homogeneous measurements of photo- 
metric parameters (effective radius, effective surface bright- 
ness, Sersic index, and luminosity of the bulge, effective ra- 
dius, effective surface brightness, concentration, and total 
luminosity of the galaxy), kinematic properties (stellar ve- 
locity dispersion and circular velocity), and masses (bulge 
mass, galaxy stellar mass, virial and dynamical mass of the 
galaxy). 

The sample demographics are as follows: 29% of the 
host ga laxies are ellipticals, 27% lenticulars, and 44% are 
spirals l|de Vaucouleurs et al.lll99ll . hereafter RC3). Regard- 
ing nuclear activity, 23% of the sample galaxies are Low- 
lonisation Nuclear Emission-line Regions (LINERs), 11% 
host H II nuclei, 25% are Seyf erts, and 8% ar e classified as 
transition objects according to lHo et al.l (|l997l ). The remain- 
ing 33% do not show central emission. 

We describe in the sub-sections below the extraction of 
all the galaxy structural parameters. 



Sky Survey (SPSS. I York et al. 



|2009| ) of the Sloan Digital 
for as many Sample A 



and B galaxies as possible. The total SDSS sample includes 
90 galaxies, 62 from Sample A and 28 from Sample B, from 
which to derive structural parameters. 

The SDSS images are already bias subtracted, flat- 
fielded and cleaned from bright stars; however, we performed 
our own sky subtraction since the SD SS pipeline sky levels 
may be fiawed for extended galaxies (jBernardi et abllioOTl : 
iLauer et aLll2007l ). This issue has been addressed in the 
SDSS -III Data Release 8 (lAihara et al.ll201ll ; lBlanton et all 
12011"). For our purposes, we estimated the sky level of each 
galaxy image by isolating five regions away from the galaxy, 
free of any contaminant, and calculating the mode of the 
sky intensities per pixel within each sky region. The average 
and standard deviation of the five sky values was then com- 
puted. The difference between our measured sky values and 
those provided by SDSS can be as large as 4%. The SDSS 
sky level is always biased high, likely due to the inclusion 
of bright, extended sources while our interactive technique 
ensures a contaminant-free selection of the sky fields. We 
find that the typical surface brightness error in the g and 
i bands is 0.1 mag arcsec"'^ at ~ 26 mag arcsec~^ and 
Hi ~ 25 mag arcsec"^, respectively (see iMcDonald et al.l 
l201ll . for more details). The large angular extent of NGC 224 
and NGC 4594 relative to the field of view thwarted their 
proper sky subtraction; these two galaxies were therefore 
excluded from our SDSS sample. NGC 221 was also dis- 
carded due to improper positioning of the image on the de- 
tector. The remaining images were fiux-calibrated based on 
the SDS S photometric zero-p oint, corrected for Galactic ex- 
tinction l|Schlegel et al.|[l998l ) a s well as for intern al extinc- 
tion and Jf-correction following IShao et all |20o3). 

The galaxy surface brightness profiles were ex- 
tracted using the isophotal fitting methods out lined in 
ICourteau etal] (fl996 ) and ,McDonald et al.l l|201ll ). These 
make use of the astronomical data reduction package 
xvistaQ. The azimuthally-averaged surface brightness 
profiles, projected onto the major axis of each galaxy, are 
shown in Figure [5] for four sample galaxies. These are the 
elliptical galaxy NGC 5127, two high- (NGC 4036) and low- 
(NGC 4477) inclination lenticular galaxies, and the spiral 
galaxy NGC 3675 which boasts a remarkably high dust con- 
tent. 

The g and i-band total magnitude of each galaxy is then 
determined by summing the fiux at each isophote and ex- 
trapolating the light profile to infinity. The g — i colour of 
each galaxy was calculated from the difference of the fully 
corrected g and i-band magnitudes. The remaining struc- 
tural parameters were measured from the i-band light pro- 
files since, of all the SDSS band passes, the i-band suffers 
least dust extinction. We extracted the isophotal radius, 
?'24.5, corresponding to the surface brightness of 24.5 mag 
arcsec^^, the half-light (or effective) radius of the galaxy, 
rc,gai, the effective surface brightness of the galaxy, ^o.gai, 
and the galaxy concentration C28 = 51og(r8o/r2o), where 
r2o and rso are the radii which enclose 20% and 80% of 



^ See jhttp : // astronomy ■ imsu . edu/holtz/xvlsta/ Index . html | 
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Figure 2. Typical g (red points) and i-band (blue points) surface brightness profiles extracted along the major axis for a few sample 
galaxies. 



the total luminosity, respectively. These i-band structural 
parameters are listed in Table [S] The g-band magnitude 
is also listed. Based on simulated models of spiral galaxies 
l|MacArthur et al.ll20o3 ). the typical error per galaxy struc- 
tural parameter is roughly 10-20%. 



3.2 Photometric Parameters of Bulges and Discs 

The structural parameters for elliptical galaxies, modelled 
typically as a single spheroid, and for spiral galaxies, mod- 
elled as the sum of a spheroid and a disc, were derived by ap- 
plying the two-dim ensional photometric deco mposition al- 
gorithm GASP2D (|Mendez-Abreu et al.|[2008l ) to the SDSS 
i— band images. 

The surface brightness of the spheroid component (typ- 



ically the entire elliptical galaxy or the bulge component for 
a dis c galaxy ) is modelled using a S ersic function (|Sersi3 
1 19681 : see also lGraham fc Drivei<l2005l ) 

7b(r) = /el0-''"«'-/-)''"-i>, (1) 

where is the effective radius, Jo is the surface brightness at 
re , and n is a shape parameter that describes the curvature 
of the radial profile. With n = 1 orn = 4, the Sersic function 
reduces to the exponential or de Vaucouleurs f unction, re- 
spectively. The coefficient 6n = 0.868 n — 0.142 ()Caon et al.l 
1993) is a normalisation term. The spheroid model ellipti- 
cal isophotes have constant position angle PAb and constant 
axial ratio Qb- 

The surface brightness distribution of t he disc compo - 
nent is assumed to follow an exponential law (|Freemanlll970l ) 
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(2) 



where h and Jo are the scale length and central surface 
brightness of the disc, respectively. The disc model elliptical 
isophotes have constant position angle PAd and constant 
axial ratio q^. The fitting algorithm GASP2D relies on a 
X minimisation of the intensities in counts, for which we 
must adopt initial trial parameters that are as close as pos- 
sible to their final values. The latter were estimated from 
our ellipse-averaged light profiles from which basic fits to 
the bulge an d disc were applied to esti mate structural pa- 
rameters fsee lMendez-Abreu et al.]|2008l . for details). The 
minimisati on is based on the robust Levenberg-Marquardt 
method bv lMore et al.l (|l98Cll ). The actua l computation ha s 
been done using the MPFIT algorithm (|Markwardtl 120091 ') 
under the IDlJj environment. 

The GASP2D software yields structural parameters for 
the bulge (7e, r-e, n, PAb, and q^) and disc (Jo, h, PAd, 
and qd) and the position of the galaxy centre (scyo)- In 
GASP2D, each image pixel intensity is weighted according 
to the variance of its total observed photon counts due to the 
contribution of both galaxy and sky, and accounting for pho- 
ton and detector read-out noise. Seeing effects were taken 
into account by convolving th e model imag e with a circular 
Moffat point spread function (|Moffatlll96i . hereafter PSF) 
with shape parameters measured from the stars in the galaxy 
image. Only the image pixels with an intensity larger than 
1.5 times the sky standard deviation were included in the 
fit. Foreground stars were masked and excluded from the fit. 
The initial guesses were adopted to initialise the non-linear 
least-squares fit to galaxy image, where the parameters were 
all allowed to vary. A model of the galaxy surface brightness 
distribution was built using the fitted parameters. It was 
convolved with the adopted circular two-dimensional Mof- 
fat PSF and subtracted from the observed image to obtain 
a residual image. In order to confirm the minimum in the 
X -space found in this first pass, two more iterations were 
performed. In these iterations, all the pixels and/or regions 
of the residual image with values greater or less than a fixed 
threshold, controlled by the user, were rejected. Those re- 
gions were masked out and the fit was repeated assuming, as 
initial trials for the free parameters, the values obtained in 
the previous iteration. These masks are useful when galaxies 
have spiral arms and dust lanes, which can affect the fitted 
parameters. We found that our algorithm converges after 
three iterations. 

The model decompositions for few elliptical galax- 
ies (NGC 4473, NGC 4636, NGC 4649) were signifi- 
cantly improved by including a dis c component, as found 
in a few other ellipt ical galaxies (|Kormendv et all l2009l : 
iMcDonald et al.]|2009l ). For the other ellipticals, the values of 
Te obtained either directly (empirically) from the light pro- 
file (in i|3.ip or by using a Sersic fitting function (from the 
two-dimensional photometric decomposition) agree within 
their respective errors. 

We eliminated 33 galaxies from our sample because of 
poor decompositio ns due to eithe r a strong central bar, a 
Freeman II profile (|Freemanlll970l ) . or just the overall inad- 
equacy of our single or double-component modelling (e.g., 

^ Interactive Data Language is distributed by ITT Visual Infor- 
mation Solutions. It is available from |http:/ /www.ittvis.com/[ 



due to the presence of strong dust lanes and/or spiral arms). 
We successfully performed photometric decompositions, as 
judged by a global x^ figure-of-merit, for the 57 galaxies ( 
38 from Sample A, 19 from Sample B) listed in Table [l] 
The latter includes the resulting bulge and disc structural 
parameters. Some examples illustrating the various fitting 
strategies are shown in Figure [3] These are the same galax- 
ies, whose azimuthally-averaged surface brightness profiles 
are shown in Figure [21 For all of them, the ellipse-averaged 
radial profiles of surface brightness, ellipticity, and position 
angle of the model image are consistent with those measured 
on the galaxy image. The differences between model and 
data found in the ellipticity and position angle of NGC 4477 
and in the ellipticity of NGC 3675 are due to the presence 
of a weak bar and strong dust lanes, respectively. These 
features are clearly seen in the galaxies' residual images. 
Nevertheless, the surface brightness residuals A/ii are re- 
markably small; |A/ii| < 0.2 mag arcsec"^ for all the galax- 
ies shown in Figure [3l except for some portions of the dust 
lanes of NGC 3675 where the residuals increase to ~ 0.4 
mag arcsec""^. Although part of the NGC 3675 disc is miss- 
ing, this does not affect the fit result. GASP2D performs a 
reliable fit as soon as the observed galaxy can be modelled 
by the sum of two axisymmetric comp onents and the field o f 
view covers at least half of the galaxy jMendez- Abreull2008l '). 

The GASP2D formal errors obtained from the min- 
imisation method are not r epresentative of the real erro rs in 
the structural parameters (|Mendez-Abreu et al.ll200^ '). In- 
stead, the estimated errors given in Table |4] were obtained 
through a series of Monte Carlo simulations. To this end, 
we generated a set of 400 images of elliptical galaxies with 
a Sersic spheroid and 400 images of disc galaxies with a 
Sersic bulge and an exponential disc; each with a different 
PA and ellipticity. The range of tested parameters for the 
simulated images was taken from the photometric analy- 
sis of nearby elliptical gala xies by i Kormendv et al.l (|2009l ) 
and face-on disc galaxies bv lGadottiT i 20091 ) (but assuming a 
wider range of axial ratios than the latter). Our tests include 
treatment for resolution effects (galaxy distance, pixel scale, 
and seeing), colour effects (accounting for V and i bands), 
inclination, and more. The two-dimensional parametric de- 
composition was applied to analyse the images of the artifi- 
cial galaxies as if they were real. The artificial and observed 
galaxies were divided in bins of 1 mag. The relative errors 
on the fitted parameters of the artificial galaxies were esti- 
mated by comparing the input and output values and were 
assumed to be normally distributed. In each magnitude bin, 
the mean and standard deviation of relative errors of arti- 
ficial galaxies were adopted as the systematic and typical 
error on the relevant parameter for the observed galaxies. 
Overall, we find that GASP2D recovers the galaxy struc- 
tural parameters with an uncertainty ranging from 1%, to 
10%, and to 20% for brighter (6.5 ^ mtot.i < 7.5), intermedi- 
ate (10.5 < mtot.i < 11.5), and fainter (12.5 ^ mtot.i < 13.5) 
sample galaxies, respectively. 

The inclination of disc galaxies was calculated from the 
fitted disc axial ratio in the i-band as 



2 ' 

% 



(3) 



wher e go is the intr i nsic d isc axial ratio. We obtain the latter 
from 'Paturel et all l|l997l) : 
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Figure 3. Two-dimensional photometric decomposition of the sample galaxies shown in Figure[2]illustrating the various fitting strategies 
adopted with GASP2D. For each galaxy, we show in the upper panels the SDSS j-band image (left), the best-fit image (middle), and the 
residual (i.e., observed-model) image (right). The lower panels show the ellipse-averaged radial profiles of the surface brightness (left), 
ellipticity (middle), and position angle (right) measured in the SDSS (dots) and model image (green continuous line). The difference 
between the ellipse-averaged radial profiles from the observed and model images are also shown. The dashed blue and dotted red lines 
represent the intrinsic surface brightness radial profiles of the bulge and disc, respectively. No disc contribution is assumed for the 
eUiptical galaxy NGC 5127. 



log go = -0.43 - 0.053 T, (4) 

where T is the galaxy Hubble type from RC3. 

The fitted i-band disc axial ratios agree well with the 
isophotal axial ratios at a surface brightness level /is = 25 
mag arcsec"^ reported in the RC3. Our disc axial ratios are 
on average 6% lower than RC3 with a standard deviation 
of 7%. We adopted the RC3 axial ratios for the lenticular 
and spiral galaxies whose GASP2D solution for the disc was 
too uncertain; this amounts to 28 galaxies in Sample A and 
another 8 in Sample B. 



was obtained from our own decomposition of the SDSS im- 
ages (Tables [3] and |4]). The maximum difference between 
our and literature values of rc, bulge is about 20%, though 
the comparison often involves different band passes which 
broadens the discrepancy. 

The aperture correction was also applied to the a mea- 
sured for NGC 2892 and NGC 5921 by Ho et al. (2009) and 
IWegner et al.l l|2003l ). respectively. 

For Sam ple B galaxies, we adopted the values of Oc 
given bv lGOgl . These were derived as the luminosity-weighted 
mean of the stellar velocity dispersion within rc. bulge. 



3.3 Stellar Velocity Dispersion 

The measured stellar velocity dispersions, o, for g alaxi es in 
Sample A were taken from the same so urces as lB0 9l. We 
applied the aperture correction of .j0rgensen et al.l (1995) to 
transform the o into the equivalent of an effective stellar dis- 
persions, (Jo, measured within a circular aperture of radius 
fe, bulge. T he ef fective radii were also taken from the same 
sources as IB09| . except for the 35 galaxies for which rc,buigc 



3.4 Circular Velocity 

The values for the galaxy circular velocity, Vc , for both ellip- 
tical and disc galaxies were collected from dif fere n t sour ces. 

We retrieved K from the compilation of I Hoi |2007|) for 
40 disc galaxies (31 from Sample A, 9 from Sample B). These 
were derived from the H I line widths available in the Hy- 
perLeda catalogue iPaturel et al.ll2003l ). The W20 and Wao 
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line widths are measured at 20% and the 5 0% of the total 
H I line profile flux. For gala:xies missing in IHoI l|2007t ) , line 
widths were found in HyperLeda for an additional 27 galax- 
ies (22 from the Sample A, 5 from Sample B). Given multiple 
sources in HyperLeda, we favoured the larger survey source 
for each galaxy in order maximise the homogeneity of the 
data basqj. All line widths were already corrected for instru- 
mental resolution. We further applied a correction for cosmo- 
logi cal stretching and bro adeni ng by gas turbulence follow - 
ing iBottinelli et all (| 19831 ) and IVerheiien fc Sancisil llooj). 
Finally, the corrected line widths W2 0.corr and W50 c orr we re 
deprojected using the prescription of iPaturel et al.! (| 19971 ) 

V20 = 0.5(10^-^*^'°^^=«— -"■'^«)/sini, (5) 
= 0.5(10^-°^^'°^^™— -"■^^'')/sini, (6) 

where the inclination, i, is listed in Table [T] and (2] We take 
the final circular velocity as the average of V20 and V50, 

K = (1^20 + V50)/2. (7) 

Following IHoI (|2007l '). we adopted a 5% error on Vc for the 
maximum velocity listed in Hyperleda. 

For 14 early-type galaxies (7 in Sample A, 7 in Sample 
B) , we adopted t he value of Vc derived from dy namical mod- 



iBinnev et al. 1990l: NGC 3115, Bender et al. 1994 


NGC 3608. [Coccato et all 20091: NGC 4314 


Quillen et al.l 


I1994I: and 9 ellioticals in iKronawitter et al.l 


2OOOI). For a 



few other galaxies (five in Sample A, one in Sample B), 
eithe r the ionised-g as (NGC 29 11, Sil'chenko fc AfanasievI 
20041: NGC 5252, iMorse et al.l Il998l). the Hi ( IC 342, 



Pizzella et all l2005l : NGC 3801. iHota et al.l [20091 ') or the 
CO kinemat ics (NGC 4526, lYoung et al.lbOOsI '). Milkv Way 



(|Baes et al.l [20(3) rotational velocity at large radii are 
assumed to represent Vc- We adopted a 10% error when the 
Vc uncertainty from models or observations was not given. 

Finally, we eliminated 7 galaxies (4 from Sample A, 3 
from Sample B) since their quoted circular velocities are 
unreahst ically low (NGC 1497, NGC 3384, and NGC 5576, 
iPaturel e t al. 2003; NGC 3642, NGC 4429, NGC 5347, and 

NGC 7052,|Hai200a)- 

The final circular velocities of 88 galaxies from Sample 
A (65) and Sample B (23) are listed in Tables [1] and [U 

3.5 Masses 

We estimated the bulge mass as Mbuigc= cx rc,buigcO"c/G, 
where G is the gravit ational constant and a — 5.0 ±0.1 
ICappellari et all 20061 ) is a dimensionless constant that de- 



pends on galaxy structure. The effective radius, bulge, 
for elliptical galaxies is extracted from their azimuthally- 
averaged light profile (Table [3]), whereas the bulge effective 
radius of lenticular and spiral galaxies is obtained from two- 
dimensional photometric decomposition (Table [4|. In deriv- 
ing Mbuige, we implicitly assumed that the measured value 
of (Te is dominated by the bulge component. For late-type 
galaxies we expect that the disc contribution to ctc results 
in an increase of the scatter of Mbuigo as a function of the 

^ We selected the most recent and highest resolution ob- 
servations available in the on-line HyperLeda catalogue 
||http : //leda ■ univ-lyoni ■ f r J up to September 2009. 



morphological type. This does not affect our analysis since 
we could measure rc, bulge (and therefore Mbuige) for only a 
few Sb-Sbc galaxies. Moreover, the ratio B /T > 0.1 except 

for three galaxies. 

iFerrarese etlo] (|2006l ) suggested a connection between 
M, and the mass of early-type galaxies calculated as 

Me,gal = a re.gaicr.VC with a = 5.0 ± 0.1. (8) 

Note that such a mass estimate is indicative of the galaxy 
mass within rc.gai and is thus an incomplete representation 
of the total galaxy mass. We calculated Me,gai for the ellip- 
tical and lenticular galaxies in our sample by adopting rc,gai 
calculated from the azimuthally-averaged light profiles (Ta- 
ble[3]|. Fo r elliptical galaxies, Mb uige = Me,gai. The virial esti- 
mator bv lCappellari et"ai] (|2006l ) captures the entire dynam- 
ical mass so long as tot al mass traces light (|Thomas et al.l 
! 20 111). I Wolf et al'l (|2010l ) found a different coeflicient (a = 4) 
for their derivation of a mass estimator for stellar systems 
supported by velocity dispersion. The latter is a good esti- 
mate of the total dynamical mass in side a radius which is 
just larger than the effective radius ijThomas et al.|[201ll ). 
The appropriate valu e of a is actually a function of the 
Sersi c shape index n ijTruiillo et al.ll2004 [ Cappellari et al.l 
|2006| ). However the exact application of the a(n) relation, 
which results in a zero-point offset for Me.g ai, does not affect 
our conclusions. We thus make use of the ICappellari et al.] 
(2006) mass estimator. Error estimates on Mbuige and Mg^gai 
account for the uncertainty on a. 

We derived the dynamical mass of the disc galaxies from 

Mdyn = RVc^/G, (9) 

where R = r24.5 (Table [3}. Most of the circular velocities 
that we derived from Hi data yield no information about 
their radial coverage. Nevertheless, Mdyn.gai corresponds to 
the galaxy mass within the optical radius, because r24.5 is 
roughly indicative of the galaxy optical radius and the obser- 
vations of spatially-resolved Hi kinematics in spirals show 
that the size of H I discs clos ely matches that of a galaxy's 
optical disc (|Ho et al.ll2008al ). 

The masses computed from Eq. ^ and Eq. @ are 
both available only for the 10 lenticular galaxies in our sam- 
ple, since they share dynamical properties of both elliptical 
and spiral galaxies. For these galaxies, (Mdyn,gai/Mc,gai) ~ 
1.27 ±1.07. 

The galaxy stellar masses, Af^^gai, were derived from 
Li,gai under the assumption of constant mass-to-light ra- 
tio {M/L)i. (M/ L)i estimates wer e inferred from our g — i 
colours following ISell et al.l l|2003l ). Their mass-to-light ra- 
tios were derived from the stellar population models of 
iBruzual fc CharlotI (2003), which are tuned to reproduce the 
ages and metallicities of the local spiral galaxies. 



4 ANALYSIS 

4.1 Correlations Between M. and Bulge and 
Galaxy Parameters 

We used the data obtained in ij3]to build scaling relations 
between M. and the bulge (i.e., the velocity dispersion, lumi- 
nosity, virial mass, Sersic index, and mean effective surface 
brightness) and galaxy (luminosity, circular velocity, stel- 
lar, virial, and dynamical mass) properties. In addition to 
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the Sample A+B, we also built a comparison sample with 
the 30 galaxies for which all the desired physical parameters 
could be measured. For each of the above parameters x, we 
assume that there exists a relation of the form 



Tl— = a + /? log — 

Mq Xq 



(10) 



where xq is a normalisation value chosen near the mean of 
the distribution of a;— values. 

The best-fit values of a and /3, their associated uncer- 
tainties, the total scatter e and intrinsic scatter eintr of the 
relation, the Spearman rank correlation coefficient p, and 
the level of significance of the correlation P were obtained 
as a function of the sample at hand. The results are given in 
Table [5l The quoted errors are all derived through a boot- 
strap technique. 

Samples A+B and the comparison sample include 
both accurate determinations and upper limits of M,. 
The proper handling of these heterogeneous data sets 
req uires that we perform a censored regress ion analy- 
sis jFeieelson fc NelsonI 1 19851: llsobe et al.l 1 19861 ) with the 
ASURVl^l package (|Lavallev et al.lll9'9i r Linear regressions 
were calculated with the EM maximum-likelihood algorithm 
implementing the techni que of expectati o n (E -step) and 
maximization (M-step) bv lDempster et al] (|l977l ) which as- 
sumes normal residuals. For each relation, we checked that 
the distribution of residuals about the fitted line was normal. 
With uncensored data, this analysis would be equivalent to 
that of a standard least-squares linear regression. ASURV 
gives as correlation coefficie nt the general ised Spearman 
rank correlation coefficient p (|Akritaslll989l) and computes 
the level of significance of the correlation P. We derived the 
scatter of the relation e as the root-mean square (rms) de- 
viation in log(M./M0) from the fitted relation assuming no 
measurement errors. This assumption does not affect our re- 
sults, since we wish to compare the relative tightness of the 
diff erent scaling relations for the comparison sample. 

IG09| focused their analysis on M, — (Jc and M,— Li, bulge- 
We have extended their analysis to a broader range of scaling 
relations by performing least-square linear regressions with 
the secure values of M, and the bulge and galaxy parame- 
ters mea sured for Sample B. W e used the MPFITExMfl al- 
gorithm (|Williams et al.ll2010l) . which accounts for measure- 
ment errors in both variables and the intrinsic scatter eintr of 
the relation. We used the IDL routine R.CORRELATE to 
compute the standard Spearman rank correlation coefficient 
p and the level of significance of the correlation. The intrinsic 
scatter of each relation was assessed by varying eintr in the 
fitting process to ensure that the reduced = 1. The total 
scatter e was derived as the rms deviation in log(M./M0) 
from the fitted relation weighted by measurement errors. 

In addition to fitting the M. -bulge and galaxy scaling 
relations as outlined above, we also look in the subsections 
below for correlations against morphological type or nuclear 
activity as listed in Tables [1] and (2] We will see that the 
latter (morphology or nuclear activity) do not play a strong 



* The FORTRAN source code of the Astronomy 
Survival Analysis Package (v. 1.3) is available at 
|http://www2. astro. psu.edu/st atcodes/asurv 

" The IDL source code of MPFITEXY is available at 



|http://purl.org/mike/mpfitexy/| 



role in any of the M.-bulge and galaxy scaling relations. 
Morphology only plays a small role in the M, — Mi^^g^i and 
in the M. — K relations. 



4-1.1 M, versus 

We have plotted the M. — cto distribution for the Sample 
A+B in Figure |4] The slope of our M. — Ue relation is 
consis tent within the err o rs wit h those bv lFerrarese fc FordI 
(|2005l ) and iLauer et al] (|2007h as discussed bv IBOGI . and 
with that bv lGOSi by default. On the other ha nd, th e inclu- 
sion of the upper limits measured by IB09| and IG09| slightly 
changes the zero point of the relation but does not appre- 
ciably affect its scatter. The M. — CTo relation is indeed the 
tightest c orrela tion that we measure; its scatter is consis- 
tent with 'G0 9| (e = . 44 ± 0.06 dex) and slightly larger 
than Terrarese k. FordI (|2005l . e = 0.34 dex). According to 
[g09| . the increased fraction of spirals may be a source of the 
increased scatter with respect to previous estimates based 
mostly on SMBH measurements in early-type galaxies. Con- 
sistently, our analysis reveals that the M. — scatter for 
late-type galaxies (e = 0.52+0.14 dex) alone is slightly larger 
than that of early-type galaxies (e = 0.38 + 0.04 dex). Never- 
theless, at small a a, some upper limits exceed the expected 
M, as the line- widths for such Iow-cto outliers are most likely 
affected by t he mass contribution of a conspicuous nuclear 
cluster (|B09r ). The scatter for the late- type galaxies reduces 
to (e = 0.38 ± 0.05 dex) by excluding the nucleated galaxies 
from the fit. Therefore, they are the cause for the observed 
difference between the scatter of early and late- type galaxies. 
We conclude that the M. — CTo relation is the same for both 
early and late-type galaxies. SO galaxies overlap ellipticals 
at high (Jo values and spirals at low Oc values. 

Barred and unbarred galaxies follow the same M, — 
Oc relation within the errors. The same is not true for 
galaxies with classical and pseudo-bulge. We have clas- 
sified our bulges according to their measured Sersic in- 
dex (n > 2 for classical an d n ^ 2 for pseudo-bulges; 
iKormendv fc Kenn icutt||2004h . We find slope differences in 
the M, — Oc relations of our 46 classical (/? = 3.86 + 0.71) 
and 11 pseudo-bulges (/3 = 5.63 + 0.86), whereas their zero 
points (a = 8.09 + 0.09 and a = 7.78 + 0.33, respectively) 
are in closer agreement within the errors. Conversely, iHrJ 
(2008) found similar slopes and different zero points. At 
face values, our and Hu et al.'s correlation coefficients are 
in agreement. This is mostly due to the large error bars on 
the M, — ae coefficients for pseudo-bulges, which translate 
into a lower level of significance of the Af, — CTc relation for 
pseudo-bulges iP < 2%) with respect to that for classical 
bu lges (P < 0.1%). This re sults supports the recent findings 
bv lKormendv et all l|201lD that there is little or no correla- 
tion between M. and pseudo-bulges. 

At high (Jo values, the M, of the Sample A+B shows a 
weak steepening with respect t o M. — Ue confirmin g pre- 
vious results bv IWvithel l|2006l ). iLauer et al] l|2007l '). and 
iDalla Bonta et all (|2o'o9h . 

Finally, Figure |3] shows the absence of trends of the 
M, — (Jo relation with nuclear activity. The different types 
of nuclear activity cover all the (Jo range, with most LINERs 
at high (Je and most Seyferts at low (Jo. This is expected 
considering the strong dependence of nuclear spectral class 
on Hubble type (and therefore on (Jo), with LINERs and 
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Figure 4. M, as a function of (Te for 94 Sample A galaxies (89 
upper limits from nebular line widths, circles; 5 upper limits from 
resolved kinematics, arrows) and 49 Sample B galaxies (squares). 
The total number of galaxies is A'^ = 143. The error bars for 
(Te are shown only in the upper panel for clarity. The lower and 
upper ends of the dotted lines correspond to M, upper limits 
estimated assuming an inclination of i = 33° and 81° for the 
unresolved Keplerian disc, respectively. Galaxies are plotted ac- 
cording to morphological type (upper pane l) and nuclear activity 
(lower panel). The dashed line is the lGOSi M, — (Te relation. The 
Sbc-|- bin includes all galaxies classified as Sbc or later. 



Seyfert nuclei being more frequent in ellipticals and spirals, 
respectively (see iHolboOSl . for a review). 



^.1.2 M, versus I/i,bu 



Igc 



The correlation between M, and the luminosity of the 
spheroidal com ponent of a galaxy is also fairly tight (see 
iGrahamI |2007| . and references therein). iGrahamI (120071 ) 

by 



M. 



b ulge 



compared the different versio n s of the 

Kormcndy fc Gc bhardj (|200 . , , 

Marconi fc Hund (j2003l 'l. and lErwin et all (|2004l ). Recently, 



McLure fc Dunlonl (l2002l ). 



Sani et al'l (|201lf ) derived the M. 



I'buige relation by mea- 
suring the 3.6 iim bulge luminosity. Since t he differing rela - 
tions do not predict the same SMBH mass, IGrahamI (20071 ) 
investigated the effects of possible biases on the M, — Lbuige 
relation including any dependency on the Hubble constant, 
the correction for dust attenuation in the bulges of disc 
galaxies, the mis-classification of lenticular galaxies as el- 
liptical galaxies and the adopted regression analysis. These 
adjustments resulted in relations which are consistent with 
each other and suitable for predicting similar M,. In par- 
ticular, a detailed photometric decomposition of the galaxy 
light profile is crucial to obtaining a representative bulge 
luminosity and therefore a reliable M, ~ Lbuige- 



Figure 5. Af, as a function of Li, bulge for 38 Sample A galaxies 
(35 upper limits from nebular line widths, circles; 3 upper lim- 
its from resolved kinematics, arrows) and 19 Sample B galaxies 
(squares) for which two-dimensional bulge-to-disc decompositions 
of the SDSS i-band images were performed. Symbols and panels 
are as in Figure |4] 



[G09| derived an updated value of the intrinsic scatter 
of the M, — Z/buige relation, which is comparable to that of 
the M, — Oe relation in ear l y-typ e galaxies, confirming the 
results of 

Im 

arconi fc Hunt! (2OO3I ). According to them, the 
scatter of the M, — Lbuige relation is significantly reduced 
when the bulge effective radius is extracted from a careful 
two-dimensional image decomposition. 

Figure [5] shows our M. — Li, bulge relation. Th e best-fit 
coefficients are in agre e ment within the errors with IGrahamI 
(.2007) an d lSani et al.1 (|201ll ). The scatter is slightly higher 
than IG09| who fit only early-type galaxies, but larger than 
that of the M. — iTc relation. All galaxies show a simi- 
lar M, — Li, bulge distribution with no dependence on mor- 
phological type or nuclear activity. We cannot test the 
claim that barred and unbarred galaxies follow different 
M. - Li. bulge relations (lGrahamll2008l: iGraham fc Lil l2009l : 
iGadotti fc KauffmannI |2009| : IHuI [2009^ ■ since the sample 
galaxies with a strong bar were excluded for photometric 
decomposition and thus lack a Li, bulge measurement. For all 
the remaining galaxies, we did not model any other compo- 
nents than the bulge and disc (see Table |4}. 



4.1.3 M, versus Mhu 



Igo 



The connection between the M, and Lbuige suggests a lin- 
ear correlat ion with mass of the spheroidal component of 
the galaxy (iMagorrian et al.lll998l: iMcLure fc Dunlodl2002l : 



iMarconi fc Huntll2003l : iHaring fc Rixll2004l ) 

The slope of our M, — Mb uige relation (Fig. |6] ) is con - 
sistent within the errors with iMcLure fc Dunlopf (|2002h . 
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Figure 6. M, as a function of M^uige for tlie same sample as in 
Figur e [5] Symbols and pan els are as in Figure |4] The dashed Uno 



is the lHaring fc Rixl ||2004| ') M.-M] 



-'bulge 



relation. 



Marconi fc Hund (|2003h . lAUer fc Richstond tOOlh and 
Sani et alj 1 20 111 ) and slightly shallower than Haring &c Rixl 



l|2004 l. This regression is not a marked improvement over 
the M, — (Te relation, as one might expect considering the ad- 
ditional fitting parameter in the mass measurement (namely 
the radius). Indeed, the M, — Mbuige is worse than the 
M, — (Tc relation. However, Mbuige is still a better proxy for 
M, than Li, bulge- Different Hubble types follow the same 
M, — Mbuige relation, with the lenticular galaxies cover- 
ing the whole range of masses. There is also no depen- 
dence of this relati o n on nuclear activity, in agreement with 
iMcLure fc Dunlod l|2002D . 



4-1-4 M, versus Sersic n and (/le, bulge) 

The two-dimensional photometric decompositions yield a 
panoply of galaxy structural parameters, including the 
Sersic shape index n and mean effective surface brightness 
(Me, bulge), both uscd as measure of the concentration of the 
bulge light. The Sersic index n has indeed been adop t ed by 
some as a good tracer for M. (|Graham et all l200ll . l2003l : 
iGraham fc Drivedl2007l ). 

Our M,—n relation (Figure O differs significantly from 
the linear and quadratic rel ations with log M. and l o g n b y 
iGraham et al] (|200ll , |2003| ) and ICraham fc DriveJ l|2007D . 
respectively. In particular, the values of n for galaxies with 
high M. are not as large as those in IGraham fc Driven 
The relation is characterised by a large scatter, small 
Spearman correlation coefficient, and low level of signifi- 
cance. Therefore, the correlation between M. and n is poor 
and the M,—n relation is not reli able for predict i ng M. . 
Our findings are in agreement with iHopkins et all l|2007bl ) 



Figure 7. M, as a function of n for the same sample as in Fig- 
ure [5] for which two-dimensional bulge-to-disc decompositions of 
the SDSS j-band images were performed. Symbols and panels 
are a s in Figure Bl The M»— n rela tions by Graham et al. (200l|, 
[loos') and 'Graham fc Drived | |2007|) are shown as the dashed and 
continuous lines, respectively. 



who found no correlation between M, and Sersic index with 
both observations and hydrodynamical simulations. 

The correlation between M, and (pe, bulge) (Figure [8| 
is also poor, lending further support to the idea that M, is 
unrelated to the light concentration of the bulge, regardless 
of galaxy morphology and nuclear activity. 



4-1-5 M, versus ii.gai 

iKormendvl l|200ll ) compared the B-band total magnitude 
with dynamically secure M. estimates in spheroids and few 
bulgeless disc galaxies. The poor correlation between M. 
and Li, gal led him to conclude that the evolution of SMBHs 
is linked to bulges rather than discs. 

We compared our own -Li,gai with M, as a function of 
the morphological type in Figure 15] The Spearman rank cor- 
relation coefficient suggest s a correlation at 17% significance 
level, which is tighter than lKormendvl (|200ll ) but not as tight 
as that for the M. — Lbuige relation. Our sample is short of 
bulgeless galaxies which may explai n the be tter agreement 
for our Al, and -Li,gai with respect to iKorme ndv (2001). On 
the other hand, total and bulge luminosity differ for disc 
galaxies. The later the morphological type, the larger the 
discrepancy resulting in both a larger slope and scatter of 
the Al, — Li, gal relation with respect to M. — Lbuigc- The 
scatter is larger and the correlation weaker when the g— band 
total luminosit y is considered. Our findings are consistent 
with iHrJ (|2009l . see their Fig. 6) and confirm that bulge lu- 
minosity is a better tracer of the mass of SMBHs than total 
light of the galaxy. 
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Figure 8. Mm as a function of (/^e, bulge) fof the same sample as 
in Figure [3 Symbols and panels are as in Figure [4] 



Figure 10. M, as a function of M^. g^i for the same sample as in 
Figure|9]for which {g — i) colours from SDSS g and i-band images 
were derived. Symbols and panels are as in Figure l4l 




10^ 10" 10'° 10" 10'^ 
Li,gai (Le) 

Figure 9. M, as a function of ii^gai for 62 Sample A galaxies 
(57 upper limits from nebular line widths, circles; 5 upper lim- 
its from resolved kinematics, arrows) and 28 Sample B galaxies 
(squares) for which azimuthally-averaged luminosity profiles from 
SDSS i-band images were extracted. Symbols and panels are as 
in Figure |4] 



4.1.6 M. versus M*,gal, Mc,gal, anrf Ai"dyn,gal 

We have tested whether the scatter in predicting M, could 
be reduced by adopting M*,gai rather than Li,gai, since 
the galaxy luminosity is a proxy for its stellar mass. The 
M, — Mj-^gai relation shown in Figure[TD]is a slight improve- 
ment over the M, — -Li,gai relation. The slopes are compa- 
rable but the Spearman rank correlation coefficient of the 
former is higher with a 3% significance level. Ellipticals and 
lenticulars follow a tighter M, — M*,gai relation than late- 
type galaxies; the disc light clearly plays an anti-correlating 
role, indicating once more that the bulge parameters drive 

SMBH correlations. 

W e adopted the mass estimator from ICappellari et al.] 
(2006) (i.e., the mass within the effective radius) for our 
ellipticals and lenticulars. The distribution of M, as a func- 
tion of Mc,gai for the early-type galaxies is plotted in Fig- 
ure 111! Th ey follow the same tr end as the M, — Mc,gai 
relation by iFerrarese et al.l (|2006l ). This is especially true 
for Me, gal > lO'^^ Mq where the fit s lope and norma l isatio n 
are consistent within the errors with IFerrarese et al.l (|2006l ). 
Less-massive galaxies seem to follow a steeper M, — Me^gai 
relation with a smaller normalisation. More data are how- 
ever needed in this mass range to carefully address any trend 
differences. Introducing a new parameter re,gai and a specific 
exponent for o"o ought to considerably reduce the scatter of 
the M, — (Tc relation but, in fact, it does not. In particu- 
lar, the scatter of the M, — relation is smaller than that 
of AI, — Mo, gal if the same sample of early-type galaxies is 
considered. 

The dynamical mass Mdyn,gai (i-e., the mass within the 
optical radius) was taken as the galaxy mass estimator for 
lenticulars and spirals. It does not correlate with M, (Fig- 
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Figure 11. A/, as a function of gai for 27 Sample A galaxies 
(24 upper limits from nebular line widths, circles; 3 upper lim- 
its from resolved kinematics, arrows) and 24 Sample B galaxies 
(squares) ranging from E to S Oa. Symbols and panels are as in Fig- 
ure |4] The dashed line is the iFerrarese et al. I 1 I2OO6I) M.-Me,gal 
relation. 



Figure 12. M, as a function of Af^yn g^i for 41 Sample A disc 
galaxies (38 upper limits from nebular line widths, circles; 3 up- 
per limits from resolved kinematics, arrows) and 6 Sample B 
(squares). Symbols and pa nels are as in Figure |4] The dashed 
line is the lHo etal] ll2008bh M.-Mdyn,gal relation. 



ure I12|l , irrespective of Hubble type and nuclear activity. 
Our ana lysis reveals that the coarse M, — Mdyn,gai relation 
found bv lHo et al.1 (|2008bl . e = 0.61 dex) does not hold when 
galaxies with Mdyn.gaiSS lO'^" M© (i.e., at the lower mass end 
of the A/dyn.gai range) are taken into account. The scatter 
of the M. -galaxy mass relation is larger and the correlation 
weaker when M*,gai and Mdyn,gai are considered. Since the 
stellar and dynamical masses include the mass contribution 
of the disc component, we conclude that the bulge mass is a 
better tracer of M, . Nevertheless, the M. — tJo is tighter and 
yet again more fundamental than the M. — Mc^gai relation. 



M. - K and K 



lations 



Assuming that K at large radii trace the dark matter 
halo, we study the possible link between SMBHs and dark 
matter halos by comparing Af. with K. A possible rela- 
tion between M, and the dark matter depends clearly on 
the radius at which K is measured. Theoretical models 
that reproduce the observed luminosity function of quasars 
llCattanc o 2001; Adams et a l. 2003; Vol ontcri ct al 20o3 : 
iHopkins et al.ll2005al : ISpringel et al.1 120051 ) suggest that M. 
scales as a power law of the virial velocity (i.e., the circu- 
lar velocity of the galactic halo at the virial radius) of the 
galactic halo of the SMBH host. However, the conversion be- 
tween Vc measured at large radii and virial velocity, depends 
on the assumed model. 

We show M, against K in Figure 1131 There is a weak 
correlation between these quantities, which is mostly driven 
by the late-type galaxies. 
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Figure 13. M, as a function of Vc for 65 Sample A disc galaxies 
(61 upper limits from nebular line widths, circles; 4 upper lim- 
its from resolved kinematics, arrows) and 23 Sample B galaxies 
(squa res) . Symb ols and panels are as in Figure|4] The dashed line 
is the|H3 ll2007l) M,-Vc relation. 
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This i s in a greeme nt with previou s studies by 
IZasov et al.l l|2005h and IHo et alj (|2008bl ). IZasov et alj 
1 20051 ) used a collection of 41 galaxies with M. , cto, and Vc 



from the literature to conclude that there is a coarse M,—Vc 
relation and that, for a given K, early- type galax ies have 
larger M, than late-type galaxies. I Ho et ahl (|2008bl ) studied 
the M, — Vc relation for a sample of 154 nearby galaxies 
comprising both early and late-type systems for which M, 
was estimated from the mass-luminosity - line width relation 
jKaspi et all I2OO0I : iGreene fc Hoi l2005l : iPeterson fc Bentj 
I2OO6I ) and Vr from H I and optical data (IHo et aLlbOOSal ). 
They found that the correlation between M, and Vc im- 
proves if only the galaxies with the most reliable Vc measure- 
ments are considered. The distribution of our spiral galaxie s 
is consistent with the M,—Vc relation bv lHo et al.l ( 2008lJ ) . 
whereas the value of M, for elliptical and lenticular galax- 
ies is almost constant over the full observed Vc range {Fig- 
ure lTHl). 

IHo et all (|2008bl ) suggested that the main source of 
scatter in the M,—Vc relation is rela ted to the dependence 
of Vc/ae on the bulge-to-disc ratio (jCourteau et al.l l2007l : 
lHdl2007l ). Sinc e the bulge-to-d isc ratio scales with the light 
concentration (iDoi et al.ill"993h . which is another way to es- 
timate the degree of bulge dominance, we derived Vc/ae as a 
function of C28. We found a general a greement with the re- 
lation found bv lCourteau et al] ([2007') . The Jeans equation 
evokes a relation between the circular velocity of galaxy and 
its velocity dispersion for a pressure -supported systera. The 
Vcjoe relation was fir st repo rted bv^ Whitmore et al.l (Il979 



and more recently by iFerr ar ose ( 200 21'). iBuvle et al 
and iPizzella et al.1 (120051 ). ICourteau et al.l (|2007l ') 

further demonstrated that the Vc—Oc relation must 
depend on a third parameter which depends on the galaxy 
structure. 

The existence of the M. — Oc relation and the absence 
of a single universal Vc — fe for all the morphological types 
is in contrast with the hypothesis that M, is more funda- 
mentally connected to halo than to bulg e, as suggested by 
iFerraresd (|2002l ') and'PizzoU a et al.l (|2005l i. We conclude that 
the SMBH m ass is assoc iated with the bulge and not the 
halo (see also |Pen^l2010l ) by analysing the 14 — CTc relation 
as a function of M, for our sample galaxies (Figure I14p . 
For a given M,, the range of Vc is on average 1.6 times 
larger than the range of ctc demonstrating that M, is driven 
by CTe and not by Vc- The few bulgeless galaxies known to 
host a SMBH (N GC 4395, iFilippenko fc IM|2003I : POX 52, 
iBarth et al.ll2004l : NGC 1042. [Shields et al.ll2008l 'l are an ex- 
ception to this scenario, which is supported by cases like 
M33. The latter is a pure disc galaxy, which doe s not show 
any evidence of a SMBH jGebhardt et al.llioOll ') but has a 



massive dark matter halo (|Corbellill200i) 



More recently, iKormendv fc Bended (|201ir ) confirmed 
our early findings (jBeifiori 201^) by extending the anal- 
ysis of the 14— Ce relation to bulgeless galaxies. They re- 
ported an absence of correlation between Vc and cTe unless 
the galaxy also contains a bulge, suggesting that the funda- 
mental driver for M, is (Tc for all Hubble types and that Vc 
plays only a small role in galaxies with a bulge. Although 
at this later cos mic time the M,—Vc i s coarser than the 
M, — (Te relation, IVolonteri et"all l|201ll ) remind us that a 
tighter connection in the past is not precluded. At earlier 
epochs, the SMBH assembly was more likely coupled to 
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Figure 14. Vc — ctc relation as a function of M, for th e same 
sample as in Figure [131 The dashed line is the IHoI l|2007l 'l Vc— o"o 
relation. 



the dark matter halo properties based on a merger-driven 
hierarchical formation scenario. Still, the scatter we mea- 
sure for the M, — Vc relation is a factor ~ 2 larger than 
that of the M, — Ue relation, which is significantly more 
than argued by IVolonteri et al.l (j201l[ ) in their analysis of 
iKormendv fc Bended (|2011^ 's data set. 



4.2 Correlations Between M, and Linear 
Combinations of Bulge and Galaxy 
Parameters 

We wish to understand whether the relations between M, 
and bulge and galaxy properties studied in H4.ll can be im- 
proved by the addition of a third parameter. For the bulge 
or the galaxy parameters x and y in Table [6l a correlation 
of the form 



M, , X n ■, V 

1°S TT- = 7 + a log — + P log — 



(11) 



is assumed, where xq and j/o are normalisation values cho- 
sen near the mean of the distribution of x and y values 
respectively. The fit parameters are the offset 7, and the 
logarithmic slopes a and p. The best-fit values of a, j3, and 
7, their associated uncertainties, the total scatter e and in- 
trinsic scatter eintr of the relation are given in Table [B] The 
Spearman rank test cannot be evaluated in the same way as 
described in H4.1l as it applies to two variables. To estimate 
the degree of correlations between the parameters in use, we 
applied the Spearman rank statistic using the projections of 
the planes once the combination of the two independent vari- 
ables has been fixed to its best-fit values. This provides an 
indication of the degree of correlation between several linear 
combinations. The Spearman rank correlation coefficient p, 
and the significance of the correlation P derived for the pro- 
jections are also given in Table [G] The quoted errors are all 
computed through a bootstrap technique. 
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For both Samples A+B and the comparison sample, 
we performed a censored regression analysis with two in- 
dependent variables assuming normal residuals. The best- 
fit parameters were computed with the EM algorithm im- 
plemented in the ASURV package which yields maximum 
likelihood estimates from a censored data set allowing, in 
the general form, N observables distributed via a multivari- 
ate Gaussian distribution. We applied this technique to the 
special case of two independent variables and another depen- 
dent variable containing the censored data. The regression 
problem solves like a least-square fit but accounting for cen- 
sored data. It uses initial guesses from an ordinary linear 
regression of the non-censored data and then converges on 
the new censored values. 

The standard deviation is estimated in the standard 
fashion, weighte d for a factor wh ich accounts for the cen- 
sored data (see llsobe et al.lll986l . for details). The gener- 
alised Spearman rank correlation coefficient, and the sig- 
nificance of the correlation for the projec tions of planes , 
are based on the technique described in lAkritasI |l98^. 
The total scatter was computed as the rms deviation in 
log(M,/M0) from the fitted relation assuming no measure- 
ment errors. 

We used our own modified version of the MPFITEXY 
algorithm to perform a least-squares linear regression of 
Sample B's data with two independent variables. Measure- 
ment errors were invoked in the derivation of best-fit pa- 
rameters, total, and intrinsic scatter. The Spearman rank 
correlation coefficient and level of significance of the cor- 
relation were also derived by using the projections of the 
best-fit plane given all possible parameter combinations. 
The actual computation was done using the IDL routine 
R.CORRELATE. 



4-2.1 M, versus Bulge Parameters 

So far we have found that M, is more strongly related 
with CTc than Li, bulge and Mbuigc. Nevertheless, the corre- 
lation between the M, and i\fbuigc suggests that a pos- 
sible combination of the bulge properties can result in a 
tighter correlation with M. . Note that the BHFP derives 
from the FP which connects re. bulge, o"e, and (M e ibuige) 
together (iDressler et all 1 19871 : iDiorgovski fc Davis! 1 19871 : 
|j0rgensen et al.l 20071 . and reference therein). Therefore, we 
have analysed the residuals of M, — a^, M, — Lbuige, and 
M, — Mbuigc relations to test whether there is the ev- 
idence for an additional dependence with an additional 
parameter. The residuals were obtained as the differ- 
ence between the observed data and their expected val- 
ues from the fitting relations in Table [5] We compared 
the residuals in log(M,/Af. [fjo]), log(A/./Af,[Li, bulge]), and 
log(M,/M, [Mbuigo]) with the bulge photometric and kine- 
matic properties (i.e., (pe, bulge), ^cbuigo, n, and (Te; Fig- 
ure lisp . There is no trend between log(M./Af. [cto]) and 
other bulge properties (Figure 1151 left column), whereas 
log(M,/Af. [Li, bulge]) and log(Af./Af. [Mbuige]) show a very 
weak dependence on To, bulge according to the the Spearman 
rank correlation coefficient and their significance intervals 
~ 6% and ~ 2%, respectively (Figure [151 midd le and left 
colum ns, respectively). This was confirmed by ISani et al] 
l|201ll ). who found no correlation of the residual with r^. 
We have tested for the more fundamental driver of Al,. 



To this aim we considered different linear combinations of 
the bulge parameters with M, (Figure [T6|. The correlation 
between A/., a^, and Tc, bulge is as tight as the M, — (Fig- 
ure [16^). Taking into account r-o, bulge does not significantly 
improve the A/, — CTe fit. The relation between M,, Li, bulge, 
and Te, bulge (Figure I16b l is not as strong as that between 
AL,, CTe, and re, bulge, but it is a slight improvement over 
M.-L i, bulge. The relation between AL,, Ce, and Li, 

bulge IS 

as tight as that between AL,, (Je, and re, bulge (Figure [TBb ). 
The relation between Af, , tJe, and Afbuige (Figure [T6hl is 
expected since Afbuige is known to correlate with AL, and it 
depends on both and re, bulge- For this reason the two cor- 
relations with AL, and either and ALbuige or re, bulge and 
Afbuige are different (but not independent) expressions of the 
same relationship. Thus, larger SMBHs are associated with 
more massive and larger bulges, as understoo d in the frame- 
work of coevolution of spheroids and SMBHs. iHopk ins et al.l 
(|2007al lbll first studied the residuals of the AL, — ALbuige and 
AL, — CTe relations with bulge properties, finding tighter cor- 
relations between the combination of the bulge parameters 
and M, . This bolstered the notion of a BHFP. The bulge pa- 
rameters Afbuige and Li, bulge are indeed physically related to 
each other through the fundamental plane relation (FP) and 
virial theorem. The existence of such BHFP has important 
implications for the largest Af, and resolves the apparent 
conflict between expected and measured values of Af, for 
the outliers in both Af, — Je and Af, — A fbuige relations (e.g. 
iBernardi eraLll2007l : ILauer et al.ll2007l ). A similar correla- 
tio n, but between AL,, re . bulge, and (/Xe, bulge), was reported 
by iBarwav fc Kembhavil |2003) for nearby ellipticals with 
measured AL, . This correlation has a smaller scatter than 
the AL, — CTe and AL, — Lbuige relations and gives further 
support to the existence of a FP-type relation for SMBHs. 

Furthermore, by comparing the tightness of the corre- 
lations between Af, and linear combination of the bulge pa- 
rameters, we argue that CTe is the parameter that drives the 
connection with AL, in the BHFP. A small contribution is 
due to re. bulge or Li, bulg e- Similar results were obtained by 
lAUer fc Richston^ l|2007l ) who compared Af, with CTe, fe, bulge 
and re, bulge- Thus, including an additional parameter to the 
relation does not improve the quality of the flt, with CTe al- 
ways being the dominant parameter. 

4.2.2 Af, versus Galaxy Parameters 

We also wish to test if the shallow Af, — Li, gai and 
AL,— AL^-.gai relations can be improved by the addition of an- 
other galaxy parameter. Therefore, we calculated the resid- 
uals of Af, — CTe, AL, — Li, gal, and AL,— ALj,,gai relations to 
look for the signature of an additional galaxy parameter. 
The residuals were obtained as the difference between the 
data points and the relation flts listed in Table [5] We com- 
pared the residuals in log(AL,/Af, [ctc]), log(AL,/Af, [Li, gai]), 
and log(Af,/Af, [Afj,,gai]) with the galaxy photometric and 
kinematic properties (i.e., C28, fe.gai, re, gai, and CTe; Fig- 
ure [T7)|. For log(AL,/AL, [CTe]), no trends with other galaxy 
properties are found whereas the log(AL,/AL, [Li, gai]) and 
log(Af,/Af, [Af*,gai]) relations show a weak dependence on re 
according to the Spearman rank correlation coefficient with 
a significance levels of 0.3% and 0.9%, respectively. We also 
examined the residuals of the AL, — (3 — i) relation to check 
that our results are not affected by the propagated errors in 
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Figure 15. Residuals of the M, — CTc, M, — Li, bulge: M^—My^^ig^ relations versus (Atc,bulgo)i ''c,bulgci '^i ""c for 38 Sample A galaxies 
(35 upper limits from nebular line widths, filled circles; 3 upper limits from resolved kinematics, open circles) and 19 Sample B galaxies 
(squares). The total number of galaxies is A'^ = 57. The Spearman rank correlation coefficient p is given in each panel. 
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Figure 16. Mm as a function of (Te and fo, bulge (a), Li, bulge and r'e, bulge (b), iTe and Li, bulge (c), Mbuige and rc, bulge (d) for the same 
sample as in Figure [151 The logarithmic slopes a and /3 and offset 7 of the fitted relation are given in each panel. 



the determination of A/^^gai. Figure [TSl shows the different 
linear combinations of two galaxy parameters with M, . A 
few tight correlations are found when considering the total 
galaxy rather than bulge parameters. The tightness always 
depends on having tJe as a fitting parameter (Table [S]), i.e., 
the Spearman rank correlation coefBcient is higher when (Te 
is taken into account. The constancy of the Spearman rank 



correlation coefficient despite the addition of other fitting 
parameters, argues that the main acting parameter is tJe. 

We conclude that the addition of bulge (e.g., re, bulge, 
ii.buige, Afbuige) or galaxy (e.g., r^gai, Li,gai, Af*,gai) struc- 
tural parameters does not appreciably modify and improve 
the M, — (Te relation. We exclude that parameter covari- 
ances would produce the same scatter in both the M, — Oc 
and M, — Oc—Tc relations. Such a conspiracy would require 
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Figure 17. Residuals of the M, — Cc, Af, — Li^gai, A/, — A/^, g^i relations versus C281 ^c,gali ^c,gali for 62 Sample A galaxies (57 upper 
limits from nebular line widths, filled circles; 5 upper limits from resolved kinematics, open circles) and 28 Sample B galaxies (squares). 
The total number of galaxies is N = 90. The Spearman rank correlation coefficient p is given in each panel. 
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Figure 18. A/, as a function of (To and r,, g^i (a), L^ gaiand ?'e,gal (b), and ii,gai (c), Af,> gai and ^e.gal (d), and Ce and Af^ gai (e) for 
the same sample as in Figure [TT] The logarithmic slopes a and /3 and offset 7 of the fitted relation are given in each panel. 



an anti-correlation between and ro, bulge (or rcgai) which 
is not observed. 



5 DISCUSSION 

The existence of the M, scaling relations implies that SMBH 
and galaxy formation processes are closely linked. Some of 
the challenges of the current models of SMBH formation 
and evolution include reproducing and maintaining the re- 



lations regardless of the sequence of galaxy evolution dur- 
ing the hierarchical mass assembly (jRobertson et al.ll2006al : 
ISchawinski et al.ll2006l : ICrotonll2009l ) . A more comprehensive 
assessment the M, scaling relations scatter thus enables a 
better characterisation of the different SMBH/galaxy for- 
mation models. 

Our large sample of galaxies with secure determination 
or upper limits of M, has enabled a thorough investigation 
of the SMBH demography over a wide range of M, , mor- 
phological type, and nuclear activity. After establishing the 
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unbiased mapping of M. upper limits against that of secure 
M, , we tested whether M, is more fundamentally driven by 
one of the several bulge (i.e., the velocity dispersion, lumi- 
nosity, virial mass, Sersic index, and mean effective surface 
brightness) and galaxy (luminosity, circular velocity, stellar, 
virial, and dynamical mass) parameters known to correlate 
with the SMBH mass, and if the known scaling relations can 
be improved with the addition of a third parameter. 

• We argue that M, is fundamentally driven by CTc, con- 
sidering that the M, — relation has the tightest scatter 
with respect to all other scaling relations. The scatter of our 
M, — (Te rel ation is comparable to I GogI but slightly larger 
than that of lFerrarese fc FordI l|2005l '). At small CTc, some M, 
upper limits exceed the expectation value as they lik ely ac - 
count for the mass of a conspicuous nuclear cluster (|B09|). 
Excluding these galaxie s from the fit reduces the scatter to 
IPerrarese fc FordI (l2005l )'s value. Barred and unbarred galax- 
ies as m iBOglandlGOgI fo llow the same M , — ae relation within 
the errors, contrary to iGrahamI (|2008l )'s findings. Classical 
and p seudo-bulges were ide ntified according to their Sersic 
index (|Kormendv fc Kennicutt 2004 ) . The M, — Oe relations 
of classical and pseudo-bulges have a different slope. The 
low significance of the M, — a"e relation for pseudo-bulges is 
in agr eement with more recent findings bv lKormendv et al.l 
l|201ll) and suggests that the formation and growth histories 
of SMBHs depend on the host bulge type. 

• The M, — -Li, bulge relation is clearl y not as tight as the 
M,—ae one. However, the fact that lCOgl found a M.—L^, bulge 
relation as tight as the M, — one is not in contradiction 
with our results, since "go? accounted only for early-type 
galaxies. The same is true for the M, — Mbuige correlation, 
although the bulge mass proved to be a better proxy of M, 
than Lbuige and it includes re. bulge as an additional fitting 
parameter. 

• C ontrary to previous fin dings (|Graham et al.l l200ll . 
I2OO3I : iGraham fc Drived 120071 ') . we find little or no cor- 
relation betw e en M, and Sersic n in agreement with 
iHopkins et all (|2007bl ). The latter stated that M, is un- 
related with the light concentration of the bulge based on 
observations and simulations. Consistently, we confirm that 
M, and (^e, bulge) are poorly correlated. 

• The correlations between M, and galaxy luminosity or 
mass are not a marked improvement over the M, — rela- 
tion. These scaling relations are strongly sensitive to the 
morphology of the host galaxies, with the presence of a 
disc playing an a nti-correlating role, as first pointed out by 
l|Kormendvll200ll '). This is a further indication that bulges 
are driving the SMBH correlations and SMBH evolution. 

• We found that the AI, — Vc relation is significantly 
coarser than that of the M, — relation, with a scatter 
about twice larger, suggesting that M, is more strongly con- 
trolled by the baryons than the dark matter. 

• To assess the need for an additional parameter 
in the M,-bulge/galax;y scaling relations, we performed 
both a residu a l analy s is and third-paramet e r fits as in 
Hopkins et al.l (l2007bl), iBarwav fc Kembhavil l|2007l '). and 
AUer fc Rich stone (20071). To this aim we considered dif- 
ferent linear combinations of bulge or galaxy parameters 
with M,. The strongest correlations always include as 
a fundamental structural parameter. The tightest relation 
is found between M,, a^, and re, bulge (also known as the 



BHFP. iHopkins et all l2007al lbh. Since its scatter does not 
change appreciably compared to that of the M, — Oc rela- 
tion, the addition of r,,, bulge is barely an improvement. This 
is a further confirmation that ctc is the fundamental param- 
eter which drives also the BHFP. 

Our findings about the tightness of scaling relations 
such as the M, — Oc may be interpr eted in the framewor k 
of self-regulating feedback in galaxies iHopk ins et al. I l|2009t ). 
These authors argued that the energy released from the ac- 
cretion of gas onto a SMBH is enough to stop further gas 
accretion, driv e away the gas, and quench star formation. 
iHopk ins et al] (|2009l ) studied the relation between M, and 
the mass of the gas accreting onto the SMBH as a func- 
tion of the distance from the central engine. The scatter 
of the M, -stellar mass relation increases at smaller radii 
(i.e., approaching the radius of influence of the SMBH) and 
decreases at larger radii (i.e., where the SMBH scaling re- 
lations are defined). This has been interpreted in terms of 
self-regulated growth of SMBHs, where the SMBH accretes 
and regulates itself without accounting for the gas supply 
but just being controlled by the bulge properties. 

In this scenario, is the property which gives the tight- 
est connection with MBH since it determines the depth of 
the loc al potential w e ll from which the gas must be ex- 
pelled. lYounger et al.l (|2008h studied self-regulated models 
of SMBHs growth in different scenarios of major mergers, 
minor mergers, and disc instabilities, to find that SMBHs 
depend on the scale at which self-regulation occurs. They 
compared the bulge binding energy and total binding en- 
ergy with M. , finding that the total binding energy is not a 
good indicator of M. mass in disc-dominated systems. This 
agrees with our flndings that the late-type systems deviate 
most signiflcantly from the Af.-galaxy and BHFP scaling 
relations. 

Several SMBH formation models predict a connec- 
tion between M. and the total mass of the galaxy 
llHaehnelt et al.lll998l:ISilk fc ReedllQOSI : lAdams et aUbOOll : 
ICroton et al.ll2006l : ICrotonll2009i rsuch that if the dark and 
baryonic matter act to form the bulge and SMBH, the dark 
halo determines the bulge and SMBH properties. There- 
fore, the mass of the SMBH and dark matter halo should 
be connected (|Cattaneoll200ll : lHopkins et al.ll2005al lbh. Nu- 
merous recent studies h ave er red along those lines. For in- 
stance, iBandara et al.l l|2009l ') studied the correlation be- 
tween M, and total (luminous-|-dark) mass of the galaxy 
estimated from numerical simulations of gravitational lens- 
ing. Their relation suggests that the more massive halos 
are more efficient at forming SMBHs than the less mas- 
sive ones; the slope of their relation suggests merger-driven, 
feedback-regulat e d pro cesses of SMBH growth. Likewise, 
iBooth fc Scha^ (I2OI0D used self-consistent simulations of 
the coevolu tion of SMBH s and galaxies to confirm the 
relation by iBandara et al.l (|2009l ) and to prove that self- 
regulat ion of the SMBH growt h occu rs on dark matter halo 
scales. IVolonteri fc NataraianI (|2009l ) investigated through 
numerical simulations the observational signature of the self- 
regulated SMBH growth by analysing the mass assembly 
history of black hole seeds. They found that the M, — Oc 
relation stems from the merging history of massive dark ha- 
los, and that its slope and scatter depend on the halo seed 
and specific SMBH self-regulation process. 
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Our observational results are not supportive of such re- 
sults. We found only a weak correlation of M. with K- This 
agrees with our expectations based on the observed tightness 
of the M. — (To relation and given the fact that the scatter of 
the Vc — (Je relation is large and morphologically-dependent 



tne Vc — Q-e relation is large and 
l|Courteau et al.ll2007l : lHo[i2007l '). 



The M,-Vc relation may be improved with homogeneous 
measurements of the K data base. Indeed, while our collec- 
tion of observed velocities is as current and reliable as possi- 
ble, especially for spiral galaxies, we still lack a fully homo- 
geneous data base. Many of us are working on improvements 
of this nefarious situation. However it is clear, as this work 
demonstrates, that the M,-a is the definitive fundamental 
relation of the SMBH-bulge connection. 
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H 


9 


33 


T 


-19 


74 
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1 
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L2 


9 


71 
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28 
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49 


2 


205.2 zt 5.8 


2 


4 
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1 


4E7 
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50 


59 


2 


— 21 


03 
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26 


2 


191.4 zt 14.4 
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2 
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s 
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NGC 3627 


SABb(s} 


T2/82 


9 
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68 
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6 
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21 
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— 20 
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2 
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T2 


12 
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10 
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66 
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9 


3 


6E7 


9 


3E6 
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46 
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69 
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65 


2 
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3 


9E8 


9 
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84 


56 


2 
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27 
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1 


1E8 
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T2 


15 


40 


2 
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71 
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63 


2 
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2E7 
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15 
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47 
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31 


2 
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81 
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60 


5 
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1 
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19 
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2 
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06 
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3 
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H 


11 


85 


2 
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71 


2 
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3 
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3 
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11 
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T2 
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79 


3 
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14 
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-19 


29 
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30 
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3 


17 
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54 


2 
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4 
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15 
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06 
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4 
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L2 


15 
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43 


5 
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14 
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Tabic 1 - Continued 
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Morph. T. 
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40 
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20.9 
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84 


88 ± 4 


28 
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SBO 
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17.0 


^09 


-19 


54 
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51 
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El 
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-16 
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56 
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Notes. — Col.(l): Galaxy name. Col. (2): Morphological type from RC3. Col. (3): Nuclear spectral class from [Hq et al.| Jl997l ). where H — H II nucleus, L — LINER, S — Scyfcrt , T — tsD 

transition object (LINER/Hli), 1 — Seyfert 1, 2 — Seyfert 2, and a fractional number between 1 and 2 denotes various intermediate types; uncertain and highly uncertain classifications '"f^ 

are followed by a single and double colon, respectively. The nuclear spectral class of galaxies marked with * is from the NASA/IPAC Extragalactic Database (NED). Col. (4): distance. 

All the distances were taken from the literature except those we obtained from V3K the weighted mean recessional velocity corrected to the reference frame of the microwave background 

radiation given in RC3. These were derived as ZJ — Vsk/^O with and Hq — 70 km s""*^ Mpc~^. For upper limits of lG09l we adopted the distances from their paper. Col. (5): References 

for Col. (4). Col. (6): Absolute corrected B m agnitude derived from (RC3) with th e adop ted distance. Col. (7): Central tj of the s tellar component within derived from the measured 

values taken from the same sources as [BOsI and corrected following |j0rgensen et al.| il995|) . The value of UGC 1841 reported by [BOsI has been corrected adopting the effective radius 

by iDonzelli et al ] 120071) . Col. (8): Incl ination of the disc galaxies derived from either the fitted disc axial ratio (from this paper) or the isophotal axial ratio at a surface brightness of ^ 
fi^ — 25 mag arcsec"^ (from RC3) and the intrinsic disc axial ratio (derived following [Paturel et al.| j 19971 ) ) . Col. (9): References for the apparent axial ratio adopted to calculate the '^'^^ 
inclination given in Col. (8). Col. (10): Circular velocity. The circular velocities were derived from either Hi line-width measurements or gas resolved kinematics available in literature O 
corrected according to the adopted inclination. Only the values marked with * are from dynamical models of the stellar component. Col. (11): References for Col. (10). Col. (12): upper 3. 

limit of lBOSi tor a Kepler ian disc model assuming i — 33° and adopted distance Col. (13): A/, upper limit of lBOSi for i = 81° and adopted dist ance. Col. (14): Af, upper limit of lGOd and 
adopted distance. 

Ref erence s. — O) [TliIIv (1988^ ), RC3, (5) Tonry et al.| <200l|) iR) iFreedman et al.| j200ll). (5 ) this paper, j^ ) [Pi zzella et al.N2005l), fvl) HyperLe da { p) iKronawitter et al .| |200C| ). @) 

[h3 j2007l "). ficll [Bhincv et al.. 1199011 . JTlh fSirdicnko AfanasievI j2004l '). jldllHota et al.l j200gl ). Ol iQuillen et al.l jl994l '). jl4ll lYoung et HI |200S|1 . JTaTlMorse et al.l jigg^ ). ^ 



Tabic 2: Properties of the galaxies of Sample B 
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4.0 


— 17.36 






158 ± 18 
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* 
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0.86 
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— 20 
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Sb 


Sl.S 
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3 
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3.4E8(2.8 ,4.9) 


NGC 3607 
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/a 
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62 zb 0.10 
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1.2E8(0.8,1.6) 
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23.0 
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— 21 


05 zb 0.10 
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57 
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2.1E8(1.4,3.2) 


NGC 3998 
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14.9 
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305 zb 15 


46 
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2.4E8(0.6 ,4.5) 
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Notes. — Col.(l): Galaxy name. Col. (2): Total extrapolated g-band magnitude of the galaxy. Col. (3): Total extrapolated i-band magnitude of the galaxy. Col. (4): Position angle of the isophotes. 
Col. (5): Ellipticity of the isophot es. Col. (6): Integ rated i-band magnitude within the isophotc at a surface brightness level of 24.5 mag arcscc"^ corrected for Galactic extinction, K-correction, 
and internal extinction following jShao et al.||2007ll . Col. (7): Radius of isophote at a surface brightness level of fj.^ — 24.5 mag arcsec"^. Col. (8): Effective surface brightness of the galaxy. Col. (9): 
Effective radius of the galaxy. Col. (10): Concentration index C28 — ^ log(^SO /'"20 ) where rgQ and r20 are the radii that enclose 80% and 20% of total light, respectively. 
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Tabic 4: Structural properties of the sample galaxies from two-dimensional decomposition of the 
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100 


76 ± 8 


75 


18.59 zb 


18 


20.66 ± 4 


08 


0.50 ± 


02 


101 


64 ± 







27 ± 


04 


11 


83 ± 0.08 


NGC 3377 


yea 






20 


41 ± 





07 


43 


52 ± 


1 


51 


3 


47 ± 





07 





70 ± 





01 


137 


05 ± 


02 






















1.00 




9 


63 ± 0.03 


NGC 3379 




bad 






































































NGC 3384 


yea 




+ cxp. 


17 


68 ± 





09 


8 


31 ± 


D 


56 


2 


33 ± 





08 





84 ± 





01 


144 


36 ± 5 


65 


19.57 zb 


11 


53.79 ± 2 


32 


0.47 ± 


01 


148 


30 ± 


09 





40 ± 


02 


10 


17 ± 0.09 


NGC 3607 


yea 






21 


14 ± 





01 


56 


34 ± 





22 


4 


70 ± 





01 





80 ± 





01 


38 


34 lb 


01 






















1.00 




9 


16 ± 0.01 


NGC 3608 


yea 






24 


20 ± 





01 


182 


19 ± 





72 


9 


03 ± 





03 





79 ± 





01 


175 


16 ± 


02 






















1.00 




9 


33 ± 0.01 


NGC 3998 






+ CX 


17 


47 ± 





13 




65 ± 





73 


2 


29 ± 





11 





84 ± 





01 


45 


48 ± 4 


74 


19.33 zb 


29 


25.51 ± 6 




0.78 ± 


03 


54 


96 ± 







45 ± 


05 


10 


80 ± 0. 14 


NGC 4026 


no^ 


bad 






































































NGC 4258 


no 


bad 






































































NGC 4261 








21 


06 ± 





01 


48 


82 ± 





19 


^ 


31 ± 





01 


^ 


77 ± 





01 


67 


96 ± 


01 






















1 00 






47 ± 01 


NGC 4342 


no 


embc 


ddcd diac 




































































NGC 4374 


yea 






20 


63 ± 





01 


63 


61 ± 





06 


4 


10 ± 


D 


01 





86 ± 





01 


37 


32 ± 


01 






















1.00 




8 


38 ± 0.01 


NGC 4459 








23 


23 ± 





01 


155 


21 ± 


D 


61 


7 


44 ± 





02 





82 ± 





01 


13 


91 ± 


01 






















1.00 




8 


78 ± 0.01 


NGC 4473 


yea 




+ exp. 


18 


05 ± 





09 


10 


60 ± 





71 


2 


23 ± 





07 





57 ± 





01 


3 


26 ± 


13 


19.64 zb 


11 


38.14 ± 1 


64 


0.56 ± 


01 


1 


46 ± 







48 ± 


02 


10 


44 ± 0.09 


NGC 4486 




light 






































































NGC 4486A 










































































NGC 4564 




embc 


ddcd diac 




































































NGC 4596 


yea 




+ cxp. 


21 


49 ± 





11 


44 


91 ± 


3 


01 


4 


43 ± 





15 





87 ± 





01 


14 


90 ± 


58 


19.48 zb 


11 


32.71 ± 1 


41 


0.67 ± 


01 


162 


85 ± 


10 





78 ± 


04 


9 


94 ± 0.09 


NGC 4649 


yea 




+ cxp. 


18 


18 ± 





09 


14 


88 ± 


1 


00 


1 


63 ± 


D 


05 





85 ± 





01 


7 


98 ± 


31 


19.03 zb 


11 


54.15 ± 2 


33 


0.78 ± 


01 


13 


78 ± 


01 





30 ± 


01 


9 


56 ± 0.08 


NGC 4697 


yea 






21 


82 ± 





01 


128 


68 ± 


D 


12 


4 


96 ± 


D 


01 





69 ± 





01 


63 


86 ± 


01 






















1.00 




8 


34 ± 0.01 


NGC 5576 


yea 






22 


42 ± 





07 


77 


58 ± 


2 


70 


8 


71 ± 





17 





69 ± 





01 


177 


21 ± 


02 






















1.00 




9 


62 ± 0.03 


NGC 5845 


yea 






17 


63 ± 





27 


4 


06 ± 





72 


3 


45 ± 





26 





66 ± 





01 


50 


79 ± 


04 






















1.00 




11 


78 ± 0.18 



Notes. — Col.(l): Galaxy name. Col. (2): Result of the two dimensional decomposition of the galaxy surface brightness distribution: yes — successful fit; no — unsuccessful fit; Col. (3): Comments 
about the fit: r^l"^ = the galaxy has been successfully fitted with a Scrsic law; r"'"''"'+exp. — the galaxy has been successfully fitted with a Sersic bulge and exponential disc; bad — inadequate 
one or two-component fit; Freeman II disc — the disc seems to follow a Freeman Type II law; embedded disc — the bulge and disc components have been fitted with a Scrsic and exponential law, 
respectively. But, the bulge is dominating the light contribution at both small and large radii; bar — a stro ng bar is prese nt in addition to the bulge and disc components: light excess — the light 1^ 
excess measured at large radii can not be parametrised with an exponential function. A cut-off profile (e.g. , [Oemled| 19761 ) is required. 

Col. (4): Effective surface brightness of the bulge. Col. (5): Effective radius of the bulge. Col. (6): Shape parameter of the bulge. Col. (7): Axial ratio of the bulge isophotes. Col. (8): Position angle of 
the bulge major-axis. Col. (9): Central surface brightness of the disc. Col. (10): Scale length of the disc. Col. (11): Axial ratio of the disc isophotes. Col. (12): Position angle of the disc major-axis. 
Col. (13): Bulge-to-total luminosity ratio. Col. (14): Total magnitude of the bulge. 
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Table 5. Fitting parameters and correlation eoeffieients foi' the different relations 







,\ 


/' 


r 














B 


49 


0.86 


< 0.1% 


8.19 ± 0.07 


4.17 ± 0.32 


0.41 ±0.06 


0.36 ±0.07 


200 km s~^ 




A+B 


143 


0.69 


< 0.1% 


7.99 ± 0.06 


4.42 ± 0.30 


0.44 ±0.05 








comparison 


30 


0.64 


< 0.1% 


7.88 ± 0.20 


4.97 ± 1.22 


0.52 ±0.10 






^i, bulge 


B 


19 


0.26 


26% 


8.60 ± 0.30 


0.47 ± 0.31 


0.56 ±0.10 


0.58 ±0.11 


10" Lq 


A+B 


57 


0.34 


1% 


8.17 ± 0.22 


0.79 ± 0.24 


0.81 ±0.13 






comparison 


30 


0.22 


24% 


7.24 ± 0.65 


0.45 ± 0.53 


0.85 ±0.32 






Afbtilge 


B 


19 


0.46 


4.6% 


8.26 ± 0.13 


0.79 ± 0.26 


0.47 ±0.06 


0.46 ±0.07 


lO" Mq 


A+B 


57 


0.62 


< 0.1% 


7.84 ± 0.12 


0.91 ± 0.16 


0.61 ±0.08 






comparison 


30 


0.44 


1.8% 


7.79 ± 0.39 


1.23 ± 0.63 


0.77 ±0.20 






n 


B 


19 


-0.1 


38% 


8.25 ± 0.18 


0.06 ± 0.80 


O.BO ±0.12 


0.61 ±0.13 


3 




A+B 


57 


0.24 


7.6% 


7.39 ± 0.19 


1.25 ± 0.60 


0.93 ±0.17 








comparison 


30 


0.30 


11% 


6.96 ± 0.36 


1.84 ± 0.89 


0.97 ±0.28 






(Me, bulge) 


B 


19 


-0.0 


46% 


8.21 ± 0.19 


1.78 ± 3.40 


0.60 ±0.10 


0.62 ±0.12 


19 mag arcsec"^ 


A+B 


57 


0.16 


23% 


7.47 ± 0.20 


4.91 ± 4.00 


0.96 ±0.16 








comparison 


30 


0.11 


53% 


6.90 ± 0.44 


7.31 ± 7.75 


1.08 ±0.35 






^gal 


B 


28 


0.44 


2% 


8.69 ± 0.22 


1.14 ± 0.40 


0.63 ±0.10 


0.62 ±0.11 


10^^ Mq 


A+B 


90 


0.34 


0.16% 


7.92 ± 0.23 


1.33 ± 0.36 


0.91 ±0.12 






comparison 


30 


0.26 


17% 


7.41 ± 0.55 


1.22 ± 0.99 


1.06 ±0.21 






*f*,gal 


B 


2S 


0.51 


0.59% 


8.21 ± 0.10 


1.43 ± 0.36 


0.49 ±0.09 


0.42 ±0.14 


lO" Mq 


A+B 


90 


0.43 


< 0.1% 


7.47 ± 0.15 


1.42 ± 0.28 


0.81 ±0.11 






comparison 


30 


0.40 


3% 


7.17 ± 0.31 


1.55 ± 0.76 


0.88 ±0.14 






Vc 


B 


23 


0.72 


< 0.1% 


7.82 ± 0.15 


3.29 ± 0.61 


0.64 ±0.08 


0.61 ±0.09 


200 km s~^ 




A+B 


88 


0.38 


< 0.1% 


6.91 ± 0.16 


4.02 ± 0.87 


0.84 ±0.11 








comparison 


30 


0.37 


4.5% 


6.82 ± 0.31 


3.59 ± 1.57 


0.84 ±0.19 






Afe.gal 


B 


24 


0.55 


0.3% 


8.17 ± 0.13 


0.83 ± 0.27 


0.54 ±0.11 


0.53 ±0.12 


10^^ Mq 


A+B 


61 


0.49 


< 0.1% 


7.86 ± 0.09 


1.00 ± 0.17 


0.49 ±0.07 




*^dyn,gal 


B 


6 


0.24 


36% 


7.62 ± 0.19 


0.32 ± 0.87 


0.46 ±0.16 


0.63 ±0.02 


10^^ Mq 


A+B 


47 


0.13 


36% 


6.10 ± 0.39 


0.55 ± 0.69 


0.96 ±0.31 







Notes. — For the variable x a correlation of the form log{i\^, / Mq) — ck + /3 log (a; /xq) is assumed for the A*" data points. For Samples A+B and the comparison 
sample, the total scatter e is defined as the root-mean square deviation in log{A/, /Mq) from the fitted relation assuming no measurement errors. For sample 
B, the total scatter takes measurement errors in both variables into account. 



Table 6. Fitting parameters and correlation coefficients for the different linear combinations 



X 




Sample 


N 


P 


P 


7 


a 


/3 


e 


eint 


XQ 


VO 


CTe 


'"e, bulge 




19 


81 


< 1% 




o.yoiEU. DO 


n ^o-l-n 01 




n ^R -l-n n7 


onn li-m a — ^ 

^uu Km S 


5 kpc 










< 0.1% 


Q 1 ^i-t-n 1 1 




n o7-i-n 1 A 

U.Z 1 iU. -1.41: 














comparison 


30 


0.68 


< 0.1% 


5.77±2.69 


4.95±1.98 


0.90±0.53 


0.37 ±0.09 








-^i, bulge 


■^e, bulge 




19 


11 


66% 


o.u.^inu.ou 


n zii -1-0 zin 


A 1 Q -1- n A A 


n fin -l-n 1 n 
u.ou inu.±u 


n R'^ -l-n 1 1 

U.OO It u. ± ± 


lOll T ^ 


J kpc 








1 0"% 


Q 1 R-l-n op; 


n TO-l-n 

U. I ^ZCU.O I 


n ni -l-n Rn 


n Qi -l-n 1 












comparison 


30 


0.20 


28% 


5.18=b4.23 


0.07+0.70 


0.71+1.14 


0.88 +0.26 








fe 


-^^i, bulge 








<, U . 1 /o 


52 QQ-l-n 1 Q 

o. oort u. ly 


Q 5252-l-n vn 

o. oolt U. / U 


n 1 i^-l-n on 
u. lorn u. zu 


n -l-n nR 
u.jy ztu.UD 


n Q53 -l-n n7 


onnb-m a — 1 


10" Lg 








"7 


f8 






'i . 'i Z zt U . ( ) o 


f 1 1 -1 -1- n '\ 


n -11 -l-n n r. 












:!() 


{\.ty.; 


' : 11. i'/ 


7.s:;±0..14 


4.!)7±2.4S 


|i.{)o±{).:;() 


O.o2 ±0.i{) 








bulge 


'"e, bulge 




19 


80 


< 0.1% 


1 -O^HJ. -LO 






n ^7 -l-n nR 


n OQ -l-n n53 


]^q11 Mq 


6 kpc 




57 


70 


< 1% 


7 1 R 

1 .oomu. ±D 


9 1 Q4-n 053 




n -l-n nA 
u . y ^ u . U4: 












comparison 


30 


0.49 


0.9% 


11.92+1.92 


2.07+0.91 


-1.46+0.83 


0.62 ±0.13 








<7e 


^bulge 










w OO-I-C"! in 
o.zzztU. iU 


o. zy ztu. o 1 


A '^o-i-n oo 

V) . o Z Zt u . z z 


n '-i7 -l-n nf; 

u . ( Zt u . uu 


n QR -l-n n7 


200 km 


10 Mq 




A + B 




70 


< () 1% 




3.83zb0.65 


U . Z / Zt U . 1 4: 


n '^o -l-n n/1 
u.oy ztu.U4i: 










comparison 


30 


0.66 


< 0.1% 


8.08+0.28 


4.09+1.84 


0.46+0.44 


0.45 ±0.10 








<7e 


'"e , gal 








< 0.1% 


Q 07-l-n 1 n 

o..^ 1 iU. J-U 




n O/i-l-n OQ 


n QR -l-n nc; 


n -l-n nR 


onn b-m — 
zuu Km s 


6 kpc 




A+B 


90 


71 


< 0.1% 


Q no-l-n no 




n nQ-l-n 01 


n An -l-n n.4 












comparison 


30 


0.64 


< 0.1% 


7.90+0.32 


4.99+1.97 


0.09+0.72 


0.52 +0.10 








-^i,gal 


'"'e , gal 








() 4% 


W 7-1 -\-<'\ 'K'i 


1 r.<:"i_|_(""i "Tir; 


— U. / ZztU . ( 4: 


( ) . f) f j zb . 1 2 


n 7 -l-n 1/1 

U. * ZtU. 14: 


10" L0 


6 kpc 




A+B 


90 


57 




o. y:Oinu. 


1 f^-l-O '^7 

o. ioinu. o 1 


— z . 1 onzu. 04: 


n 77 -l-n in 

U. 1 1 ZtU. lU 










comparison 


30 


0.58 


0.2% 


7.73+0.75 


3.06+1.79 


-3.69+2.03 


0.97 ±0.20 








<7e 


^i,gal 








< 1% 


53 07-t-n on 


/I no-t-n SI 


U. iZit u. * 


n "^7 -l-n ni^ 

U. i zt U. UO 


n -l-n n 
U.oO ^u.u 


onn L-m c — 1 
zuu Km. s 


10" Lq 




A+B 


90 


70 




u n/i -i-fi 1 

O.U41:ltU. io 


A 7 /I -l-n R 1 


A ni"i-i-n O'^ 
u. Uoztu. Zo 


n /1 1 -l-n n/i 

U.'ii ztU.U'l 












comparison 


30 


0.64 


< 0.1% 


7.90+0.36 


4.95+1.99 


0.06+0.52 


0.52 ±0.11 










^e,gal 


B 


28 


0.62 


< 0.1% 


7.53+0.27 


3.85+0.84 


-3.19+1.03 


0.68 ±0.16 


0.58 ±0.15 


10^^ Mq 


5 kpc 




A+B 


90 


0.59 


< 0.1% 


7.29+0.15 


2.74+0.41 


-2.33+0.51 


0.66 ±0.09 










comparison 


30 


0.63 


< 0.1% 


6.74+0.51 


3.07+1.65 


-3.51 + 1.82 


0.78 ±0.19 








CTe 




B 


28 


0.83 


< 0.1% 


8.21 + 0.08 


3.59+1.23 


0.32 + 0.55 


0.37 ±0.05 


0.34 +0.05 


200 km s~^ 


10^^ M0 






A+B 


90 


0.70 


< 0.1% 


8.02+0.07 


4.68+0.56 


0.08+0.24 


0.41 ±0.58 










comparison 


30 


0.68 


< 0.1% 


7.93+0.21 


4.28+1.79 


0.61+0.54 


0.48 ±0.11 









Notes. — For the variable (x,y) a correlation of the form log(M,/MQ) = 7 + ck log(a[;/cco) + ^log(y/yo) assumed for the AT data points. For Samples 
A+B and the comparison sample, the total scatter e is defined as the root-mean square deviation in log(M, /Mq) from the fitted relation assuming zero 
measurement errors. For sample B, the total scatter takes measurement errors for all the variables into account. 
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